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Abstract 



We give a physical description in terms of percolation theory of the phase 
transition that occurs when the disorder increases in the random antiferro- 
magnetic spin-1 chain between a gapless phase with topological order and a 
random singlet phase. We study the statistical properties of the percolation 
clusters by numerical simulations, and we compute exact exponents charac- 
terizing the transition by a real-space renormalization group calculation. 
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Antiferromagnetic quantum spin chains are known to present vastly different physical 
properties according to whether the spin value is integer or half-integer Q. In the case of 
integer spin, there is a (Haldane) gap for magnetic excitations and short-range magnetic 
order. In the case of half-integer spin, there is quasi-long range spin order with algebraically 
decaying spin correlations and no gap. An appealing physical picture of the integer-spin case 
is given by the so-called valence-bond-solid (VBS) wavefunctions 0. With these approxi- 
mate ground states, it is easy to see that there is a hidden long-range order that is non-local 
in terms of the true spins-1 [[3]. These properties are well established for the pure systems 
but it is obviously important to understand what happens in the disordered case. Bond 
disorder in a quantum spin chain may be realized experimentally by chemical substitution 
of the ligand that induces superexchange between the spins, for example. This leads to a 
Heisenberg Hamiltonian: 

H — Si ■ S i+1 , (1) 

i 

where Si are quantum spin operators and the couplings J, are quenched random variables. 
In the spin-1/2 case, it is known 0|| that for any strength of the disorder the system goes 
to the so-called random singlet phase M in which the spins are paired into singlets over 
arbitrarily large distances. The effect of randomness on other quantum spin chains has also 
been investigated recently |7|-pl|. 



In the spin-1 case, it has been shown [ 12 , 1 3 1 that this random singlet phase appears only 



when the randomness is strong enough. These results have been obtained by an asymp- 
totically exact real-space renormalization group (RG) study on an effective spin-1/2 model. 
For weak disorder, there is a phase which is gapless but sustains hidden long-range order. 
This is analogous to gapless superconductivity that occurs when doping a superconductor 
with magnetic impurities; here the role of the superconducting condensate is played by the 
hidden order. 

In this Letter, we use a real-space renormalization procedure that leads to an appealing 
percolation picture of the phase transition that occurs in the disordered antiferromagnetic 
spin-1 chain. We formulate the hidden long-range order in terms of a macroscopic VBS- 
cluster. This allows us to define several critical exponents to describe this phase transition. 
We study them by direct numerical renormalization and compare to the exact values that 



we obtain for the effective model of Ref. [12 which is in the same universality class. 



An efficient technique to study random spin chains is the real-space renormalization 
group @||. The basic idea is to integrate out the strongest bond of the chain and to 
compute in perturbation theory the corresponding renormalization of the coupling strengths 
distribution. This can be done iteratively and, in the case of the S=l/2 chain, leads to a 
universal fixed point distribution which is extremely broad, the so-called random singlet 
phase Q. However, when generalized to higher spin values, this method is blind to the 
difference between integer and half-integer spins. To capture the physics of the Haldane 
phase of the random spin-1 chain, we need to introduce an extended set of spin chains 
described by the Hamiltonian (p]), where Si is a spin-operator of size Si = 1/2 or Sj = 1, and 
where the couplings { Jj} can be either ferromagnetic (F) or antiferromagnetic (F), but have 
to satisfy the following constraint : for any chain segment {i,j}, the classical magnetization 
must satisfy \m it j\ < 1. 
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This condition for j = i + 1 implies that there are four types of bonds : i) F bond between 
two spins S=l/2, ii) AF bond between two spins S=l/2, hi) AF bond between one spin S=l 
and one spin S=l/2, iv) AF bond between two spins S=l. The decimation procedure is 
the following : to each bond (z), we associate the energy difference A, between the highest 
state and the lowest state of the Hamiltonian JiSi ■ Si + \. We pick up the bond (Si, Si+x, J^j 
corresponding to the strongest Aj of the chain. To define the renormalization rule for this 
bond, we divide the four-spin Hamiltonian into Hi = ho + hi where ho = JiSi ■ Si+i and 
hi = Ji^iSi-i ■ Si + Ji + \Si + \ ■ Si + 2- We treat h\ as a perturbation of ho to find the effective 
Hamiltonian replacing Hi when the highest energy state of ho is removed. This leads to the 
four following renormalizations rules: 
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These rules are a slight modification of those proposed by Hyman |I3|| . Our renormal- 
ization procedure is entirely consistent from the point of view of the progressive elimination 
of the highest energy degrees of freedom : all the energy scales of the new bonds are always 
smaller than the energy scale of the old bond. It is also easy to check that this procedure 
is "closed" inside the enlarged set of spin-chains defined above. In particular, spins higher 
than 1 cannot appear in this scheme [14 . 

If we represent each spin-1 as the symmetrization of two spin-1/2, the rules 2, 3, and 4 
for AF bonds may then be interpreted as the formation of a singlet between two constitutive 
spin-1/2. The rule 1 for a F bond between two spin-1/2 corresponds to their symmetrization. 
As a consequence, at the end of the renormalization procedure, when there are no free spins 
left, the chain decomposes itself into a set of disconnected clusters that have the structure of 
a VBS-state (see Figure 1). The topological order can be probed by use of the string order 
parameter defined as: 
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FIG. 1. The ground state of the random chain obtained from the real-space renormalization. 
Each spin S=l is written as two symmetrized spins S=l/2: the black dots related by a vertical 
thick line. Spins are grouped in clusters related by singlet bonds that are the thin lines. There are 
two clusters of size 2, one of size 4, and one larger cluster that extends beyond the limits of the 
figure. 



For the pure VBS state on a finite ring of n spins S=l, one has t^- = 4/9 (up to corrections of 
0(3~ n ) that we neglect in the following). Now for the disordered system, one has ty =4/9 
if the two sites do belong to the same cluster and Uj = if they don't. We thus consider, 
for a given realization of the disorder on a chain of N spins, the spatial average of the string 
order parameter Y,i,jtij/N 2 which is equal to 4/9 x T where: 

c i,j 

Here c is the cluster index and n c is the number of spins of the cluster, i.e. its size (which 
is not directly related to its spatial extent). This quantity T is the probability that two 
randomly chosen spins belong to the same cluster. It is also the mean size of the cluster 
containing a randomly chosen spin, divided by N. This order parameter can be non-zero in 
the large-iV limit only if at least one VBS-cluster contains a finite fraction of the spins of 
the chain. 

We have simulated the renormalization procedure on spin-1 chains containing iV sites 
with periodic boundary conditions. The initial couplings Jj are uniformly distributed in the 
interval [1, 1 + d]. The parameter d measures the strength of the initial disorder. We have 
numerically implemented the renormalization rules on M initial independent samples, and 
averaged quantities over these different realizations of the disorder. Typically the averaging 
process has been performed on N x M « 10 9 spins. The plot of T as a function of the 
disorder for various sizes in Figure 2 shows clearly that the phenomenology for the clusters 
is similar to the physics of percolation. Indeed we observe a small-disorder phase in which 
there is exactly one macroscopic cluster (c=l) that contains a nonzero fraction n\/N of the 
spins in the thermodynamic limit. There is also a distribution of finite clusters, i.e. non- 
diverging with N. This phase with T ^ has long-range hidden order. There is a critical 
value d c of the disorder for which the diverging cluster disappears, i.e. T vanishes. Then, 
for larger randomness d > d c , there are only finite clusters and T = 0. 

Following the notations of percolation theory WE\ , we denote by (3 the exponent describing 



the vanishing of the fraction ni/N of spins in the macroscopic cluster. The order parameter 
thus scales as T ~ (d c — d) 2f3 for d < d c . The finite size scaling study of Fig. 2 leads to the 
following estimate: 

2/? = 1.0(1). (4) 
On Figure 3, we have plotted the mean size of finite clusters, which plays the role of the 



4 




disorder d 

FIG. 2. String order parameter T as a function of the disorder d for sizes N = 2 14 — 2 20 . The 
critical point is d c = 5.76(2) 
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FIG. 4. Cluster size distribution m c {s) as a function of size s, at criticality. The slope leads to 
r = 2.2(1). 



percolation susceptibility: 



= 

X ~ L> N' 

c>l ly 



(5) 



In the sum, note that we omit the biggest cluster (c=l). Denoting by 7 the exponent 
describing the divergence of this quantity, \ ~ \d c — d|~ 7 , we have performed a finite size 
scaling study which leads to the value : 
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1.2(1). 



(6) 



Let us call rrid(s) the number of clusters of size s, divided by N. Then standard scaling 
leads to the formula valid in the limit N 



00 p| 



m d {s) ~ — / \{d c - d)s c 
s T 



(7) 



where f(z) is a scaling function. The exponent a defines the divergence of the characteristic 
cluster size as \d c — d\~ l / a and the exponent r characterizes the algebraic decay of the 
distribution of cluster sizes at the critical disorder d c . These exponents can be related to (3 
and 7 as in percolation theory according to (3 — (r — 2)/a and 7 = (3 — r)/a. In Fig. 4, we 
have plotted the distribution of cluster sizes at criticality d = d c . Measurement of the slope 
leads to the following estimate of the exponent r: 



r = 2.2(1). 



(8) 
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We now show that the exponents that appear naturally in the percolation picture can be 
obtained exactly. First, we explain how the VBS-clusters appear in the random dimerized 
spin-1/2 chain model of Hyman and Yang [12], which is an effective model for the random 



AF spin-1 chain. We briefly recall the properties of the effective model that are necessary 
for our purposes. In this model, all even bonds are AF, whereas the odd bonds are either F 
or AF. At any stage of their renormalization procedure, the energy scale fl is given by the 
strongest AF bond of the system, so that the odd bonds separate into two groups : group 
A contains all AF bonds and F bonds weaker than Q, while group B contains all F bonds 
stronger than Q. In this effective model, perfect VBS order means singlets over all even 
bonds of the initial chain. Indeed, at the stable fixed point of the Haldane phase found in 
HI by solving the renormalization flow equations, all odd bonds are much weaker than even 



bonds so that only singlets over even bonds are generated by decimation. It is convenient 
to introduce in their procedure an auxiliary variable fx for each odd bond still surviving at 
scale Q. The variable fx is by definition the number of singlets already made over even bonds 
of the initial chain that are contained in this odd bond. The evolution rules for this new 
variable are the following : when an odd bond of variable (fx) is decimated, a finite cluster 
of size (fx + 1) is terminated; when an even bond surrounded by two odd bonds of group 
B of variables fi\ and fx 2 respectively, a finite cluster of size (fx\ + fx 2 + 2) is terminated; 
when an even bond surrounded by two odd bonds of group A, or surrounded by one odd 
bond of group A or one odd bond of group B, with respective variables fi\ and fx 2 , the new 
odd bond generated by this decimation inherits the variable (fi = fx\ + fx-i + 1). Following 
the method outlined by Fisher ||, we find that, at the fixed point describing the transition 
of the effective model, the auxiliary variable fx scales with T = In (Qq/Q) (where Qo is the 
initial cut-off) as : 

fx oc with <p = V5. (9) 

This has to be compared with the scaling I oc T 3 of the auxiliary variable / that counts the 
number of initial bonds in a surviving bond at scale T. Deviations from the critical point 



are driven by a relevant perturbation [12] that scales as T + with A + = (v 13 — l)/2. As a 



consequence, the exact exponent for the string topological order parameter (gD is |L6[ : 

w= 2(3-0 = 4(i <.- y/5) = 1 17278 



A+ V13- 1 

the exponent of the percolation susceptibility is : 

(2^-3) = 2(2^5-3) = 

and the exponent r of the scaling form (^) for the distribution of cluster sizes : 

3 3 

r = 1 + - = 1 + -= = 2.34164..., (12) 
in agreement with our numerical estimates (f|), (0), (H). 
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In this Letter, we have given given a percolation picture of the random spin-1 chain. 
Topological long-range order characteristic of the Haldane phase is due to the existence of 
a single cluster of spins related in a VBS manner whose size diverges in the thermodynamic 
limit. The transition to the random singlet phase with increasing disorder leads to the 
vanishing of the string order parameter with a exponent 2(3 = 1.17 while the percolation 
susceptibility diverges with an exponent 7 = 1.13. These exponents are exact because the 
RG procedure is asymptotically exact pfl. 
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